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A  command  procedure  has  been  developed  for  the  U.S.  Army 
Research  Laboratory  (ARL)  molecular  beam  epitaxy  (MBE)  computer 
control  system  that  allows  a  user  to  set  up  the  system  for  an 
automated,  unattended  start  each  morning.  The  automated  sequence 
consists  of- 

1.  A  system  safety  check  to  determine  if  cell  ramping  should  be 
allowed. 

2.  A  cell  temperature  ramp  to  an  outgassing  temperature. 

3.  An  outgassing  of  cells. 

4.  A  ramp-down  of  cells  to  nominal  operating  temperatures. 

5.  An  automated  setup  through  an  iterative  process  of  flux 
measurements  and  changes  of  temperatures  until  desired  targets 
are  reached. 

This  command  procedure  allows  a  daily,  safe  start-up  of  the  MBE 
system  and  generates  identical  flux  settings  that  improve  the  crystal 
growth  reproducibility.  Typically,  one  can  save  two  hours  or  more  of 
a  work  day  by  using  this  automated  procedure. 
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1.  Background 


The  U.S.  Army  Research  Laboratory  (ARL)  molecular  beam  epitaxy  (MBE) 
system  is  controlled  by  a  PC-based  system  called  "Molly,"  which  is  sup¬ 
plied  by  EPI  Technologies,  Inc.  Molly  provides  a  script  language  that  can  be 
used  to  create  command  procedures  that  execute  customized  sequences  of 
actions  on  the  MBE  system.  Possible  actions  are  reading  and  setting  cell 
temperatures,  opening  and  closing  shutters,  reading  pressure  gauges,  and 
turning  the  substrate  holder. 

The  PC-based  control  system  replaced  an  older  PDP-11  system.  The  short¬ 
comings  of  this  latter  system  were  clarified  after  a  malfunction  during  which 
the  MBE  machine  had  been  programmed  to  start  a  cell  up-ramp  when  the 
liquid  nitrogen  had  been  inadvertently  turned  off.  This  caused  damage  to 
the  growth  system,  requiring  venting  and  replacement  of  some  cells.  Al¬ 
though  the  PDP-ll-based  system  could  read  pressures,  it  did  not  allow  de¬ 
cisions  to  be  programmed  in  to  make  actions  conditional  on  any  system 
status  parameter.  Molly  provided  a  solution  to  this  problem  but  required 
custom  written  code.  This  report  describes  the  result  of  that  effort. 
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2.  Flow  Chart 


The  logic  behind  the  developed  script  is  that  the  MBE  system  is  idling  over 
night,  with  the  evaporation  cells  at  a  low  temperature  at  which  the  evapo¬ 
ration  rate  is  negligible.  After  the  system  vacuum  is  checked  to  ensure  a 
safe  up-ramp,  the  cells  are  slowly  brought  up  to  an  outgassing  temperature 
above  the  estimated  set  points  for  growth.  The  up-ramp  is  typically  0.5  hr, 
to  allow  the  cells  to  thermalize  to  avoid  stresses.  At  the  peak  temperature, 
the  shutters  are  opened  for  about  10  min.  to  allow  material  that  may  have 
condensed,  at  or  near  each  cell  at  the  idling  temperatures,  to  be  evaporated 
so  as  to  provide  cleaner  molecular  beams  during  growth.  After  outgassing, 
the  shutters  are  closed  and  the  cell  temperatures  are  lowered  to  the  previ¬ 
ous  day's  set  points. 

Because  material  is  consumed  during  growth  and  the  temperature  sensor 
in  each  cell  does  not  perfectly  represent  the  melt  temperature  from  day  to 
day,  the  previous  day's  set  points  typically  do  not  exactly  reproduce  the 
previous  day's  fluxes.  The  set  points  must  therefore  be  changed  based  on 
the  difference  between  measured  fluxes  and  the  target.  All  cells  obey  a  lin¬ 
ear  relationship  between  the  logarithm  of  the  flux  and  the  inverse  of  the 
absolute  temperature  of  the  cell.  This  relationship  is  used  to  calculate  the 
needed  temperature  change  based  on  the  measured  flux  difference.  After  a 
new  temperature  is  set,  the  computer  is  programmed  to  wait  a  predeter¬ 
mined  time  to  let  the  cell  reach  equilibrium  before  a  new  measurement  is 
taken.  Some  hysteresis  is  typically  experienced  in  this  process  that  requires 
up  to  about  six  repetitions  before  acceptable  accuracy  is  reached.  The  accu¬ 
racy  ( I  (target-measured) /target  I )  is  a  variable  that  is  typically  set  to  0.0025 — 
a  level  of  precision  for  which  a  human  operator  seldom  can  muster  the 
patience. 

The  actual  process  of  measuring  the  fluxes  has  been  designed  to  avoid  the 
flux  transients  typically  seen  when  shutters  are  opened.  These  transients 
are  caused  by  the  fact  that  with  the  shutter  in  closed  position,  heat  from  the 
melt  surface  is  radiated  back  from  the  shutter  to  the  melt.  When  the  shutter 
is  abruptly  opened,  the  steady  state  is  interrupted  and  heat  radiation  is  lost 
from  the  melt  at  a  higher  rate,  resulting  in  a  drop  in  the  flux.  After  some 
time,  the  thermocouple  at  the  bottom  of  the  crucible  experiences  a  drop  in 
the  melt  temperate,  prompting  the  controller  to  increase  the  power  to  the 
cell  until  the  temperature  set  point  is  restored.  After  this  control  sequence 
has  reached  a  new  steady  state,  the  flux  is  stable.  The  measurement  of  the 
flux  must  consequently  be  done  at  this  point  or  later  and  not  during  the 
transient. 

The  ionization  gauge  used  for  flux  measurements  is  turned  away  from  the 
cells  when  it  is  not  used  for  measurements  to  increase  its  lifetime.  A  flux 
measurement  sequence  thus  consists  of  (1)  opening  the  shutter  for  a  prede¬ 
termined  time  (usually  10  min.),  (2)  turning  the  gauge  toward  the  cell,  (3)  av¬ 
eraging  of  10  flux  readings  to  determine  the  flux  plus  the  background  pres- 
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sure,  (4)  closing  the  shutter,  (5)  waiting  for  the  gauge  reading  to  stabilize, 
(6)  averaging  of  10  flux  readings  to  determine  the  chamber  background  pres¬ 
sure,  (7)  subtracting  the  second  reading  from  the  first  to  obtain  the  net  flux, 
and  (8)  turning  away  the  gauge  from  the  cells  again.  If  the  measured  and 
target  fluxes  deviate  more  than  the  preset  accuracy,  a  new  temperature  is 
calculated  and  set.  The  system  then  waits  for  the  cell  to  stabilize  at  the  new 
temperature. 

During  the  calibration  sequence,  the  pressures  and  temperatures  are  logged 
in  a  standard  file.  The  data  in  this  file  can  be  extracted  and  plotted  as  in 
figure  1  on  page  4.  A  good  practice  is  to  save  the  log  file  with  the  day's  date 
for  future  reference. 

For  the  script  file  to  properly  execute,  certain  information  must  be  provided. 
The  file  (Calib.cmd)  can  be  opened  with  any  text  editor.  I  have  used  WinEdit, 
which  is  a  shareware  program  editor.  The  script  file  has  been  written  with 
enough  comments  next  to  the  variable  declarations  to  indicate  what  needs 
to  be  entered.  Typical  inputs  are — 

•  Which  cells  are  to  be  ramped. 

•  To  what  temperature  the  cells  are  to  be  ramped. 

•  If  the  cells  are  or  are  not  to  be  included  in  the  flux  calibration  sequence. 

•  What  the  flux  targets  are. 

•  The  start  time  of  the  ramps. 

Optionally,  other  parameters  can  be  changed,  although  the  default  values 
normally  provide  good  performance.  These  parameters  include  the  calibra¬ 
tion  accuracy,  the  slope  and  intercept  of  the  flux  versus  temperature  lines, 
the  outgassing  temperature  expressed  as  a  percentage  above  the  growth  set 
point,  and  the  chamber  pressure  that  cannot  be  exceeded  if  ramping  is  to  be 
started.  The  Appendix  contains  the  script  listing. 
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3.  Summary  and  Conclusion 


Control  code  has  been  developed  to  allow  unattended  start-up  of  an  MBE 
system.  The  code  has  been  used  and  tested  numerous  times  and  delivered 
very  accurate  growth  parameters  and  subsequent  crystal  layers  with  ex¬ 
tremely  small  thickness  and  composition  errors.  In  at  least  one  case,  the  cell 
up-ramp  was  not  started  because  of  an  excessively  high-growth  chamber 
pressure,  thus  preventing  potential  harm  to  the  system. 


5 


Appendix.  Script  Listings 


The  following  is  a  script  listing  of  commands  that  execute  a  customized 
sequence  of  actions  on  the  MBE  system. 

/* 

/♦★★★♦♦★♦a******************************************************************************/ 


/*  */ 

/*  This  command  file  performs  a  growth  chamber  pressure  check,  a  cell  up-ramp,  */ 

/*  a  cell  degas,  and  flux  calibration  starting  at  a  predetermined  time.  */ 

/*  */ 

/*  Stefan  Svensson,  ARL  May  7,  1997  */ 


/***************************************************************************************/ 

/* 

#include  <stdlib.h> 
flinclude  <signal.h> 

# include  <unistd.h> 

#include  <cells.h> 

#include  <shutters.h> 

# include  <mbe.h> 

/**/ 

/*  ======-===  Declarations  and  initial  values  ====“=-==*/ 

/**/ 


/*  START  BY  FLAGGING  THE  CELLS  TO  BE  RAMPED  */ 


int 

rampGa 

= 

1 

/* 

Ramp 

flag 

for 

Ga 

(l=yes 

0 

=  no) 

*/ 

int 

rampAl3 

= 

0 

/* 

Ramp 

flag 

for 

A13 

(l=yes 

0 

=  no) 

*/ 

int 

rampAl4 

= 

1 

/* 

Ramp 

flag 

for 

A14 

(l=yes 

0 

=  no) 

*/ 

int 

rampln 

= 

1 

/* 

Ramp 

flag 

for 

In 

{ l=yes 

0 

=  no) 

*/ 

int 

ramps i 

= 

1 

/* 

Ramp 

flag 

for 

Si 

(l=yes 

0 

=  no) 

*/ 

int 

rampBe 

= 

0 

/* 

Ramp 

flag 

for 

Be 

(l=yes 

0 

=  no) 

*/ 

int 

rampSb 

= 

0 

/* 

Ramp 

flag 

for 

Sb 

(1-yes 

0 

=  no) 

*/ 

/*  THEN  ENTER  THEIR  SET  POINTS  (they  will  be  outgassed  at  a  2.5%  higher  temp)  */ 


double  setpGa  = 

956.3 

/*  Target  temp  for  Ga 

ADJUST 

*/ 

double  setpAl3  - 

1135.2 

/*  Target  temp  for  A13 

BASED 

*/ 

double  setpAl4  = 

1111.1 

/*  Target  temp  for  A14 

ON 

*/ 

double  setpln  = 

790.7 

/*  Target  temp  for  In 

PREVIOUS 

*/ 

double  setpSi  = 

1327.4 

/*  Target  temp  for  Si 

CAL 

*/ 

double  setpBe  = 

922.5 

/*  Target  temp  for  Be 

*/ 

double  setpSb  = 

400.0 

/*  Target  temp  for  Sb 

*/ 

/  *  THEN  FLAG 

THE  CELLS  WHICH  WILL  BE  FLUX  CALIBRATED 

*/ 

int  calGa  -  1; 

/*  Calib  flag  for  Ga 

( l=yes 

0  = 

no) 

*/ 

int  calAl3  =  0; 

/*  Calib  flag  for  A13 

( 1  =yes 

o  = 

no) 

*/ 

int  calAl4  =  1; 

/*  Calib  flag  for  A14 

( l=yes 

0  = 

no) 

*/ 

int  cal In  =  1; 

/*  Calib  flag  for  In 

( l=yes 

0  = 

no) 

*/ 

int  calSb  =  0; 

/*  Calib  flag  for  Sb 

<l=yes 

0  = 

no) 

*/ 

/*  THEN  ENTER  THE  FLUX  TARGET  VALUES 

*/ 

double  targetGa 

=  5.25E-7;  /*  Target  flux  for  Ga 

*/ 

double  targetAl3 

—  1.062E-7;  /*  Target  flux  for  A13 

*/ 

double  targetA!4 

-  0.844E-7;  /*  Target  flux  for  A14 

* 

double  targetln 

=  1.63E-7;  /*  Target  flux  for  In 

* 

double  targetSb 

=  1.0E 

-7;  /*  Target  flux  for  Sb 

* 

/*  AND  THE  PRECISION  OF  THE  CALIBRATION  (cal  ends  when  abs (  ( flux- target ) /target  )  <  prec  */ 


double  precGa  =  0.0025; 
double  precA13  =  0.0025; 
double  precAl4  =  0.0025; 
--double  precln  =  0.0025; 
double  precSb  =  0.0025; 


/*  Target  presision  for  Ga  */ 

/*  Target  presision  for  A13  */ 

/*  Target  presision  for  Al4  */ 

/*  Target  presision  for  In  */ 

/*  Target  presision  for  Sb  */ 
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Appendix 


/* 


FINALLY,  DECIDE  WHEN  TO  START  THE  EXECUTION 


/*  Use  the  formula:  */ 
/*  Finish  time  =  Start  time  +  45  +  30*number  of  cells  to  cal  (min)  */ 
/*  (The  actual  time  depends  on  accuracy  of  initial  set  point  and  desired  precision  */ 


int  hour  =  8 ; 
int  minute  =  32  ; 


/*  Start  hour  use  24  hr  clock  (0-23  valid)  */ 
/*  Start  minute  */ 


DON'T  FORGET  TO  PUT  SHUTTERS  AND  CAR  IN  REMOTE 


END  OF  STANDARD  ENTRIES 


double  backtest  =  2.E-9; 
double  ramptime  =  30.0; 
double  transient  =  150.0; 
double  outgastime  =  10.0; 
double  outgas  =  2.5; 

double  slopeGa  =  -12670.0 
double  slopeAl3  =  -10000.0 
double  slopeAl4  =  -12670.0 
double  slopeln  =  -15707.0 
double  slopeSb  =  -10000.0 

double  stabGa  =  900.; 
double  stabAl3  -  600.; 
double  stabA!4  =  600.; 
double  stabln  =  600.; 
double  stabSb  =  600.; 


/*  Background  pressure  for  up-ramp  test 
/*  Length  of  cell  ramp  (min) 

/*  Wait  time  to  avoid  shutter  transient 
/*  Outgas  time  after  up- ramp  (min) 

/*  Percent  temp  increase  to  outgas  at 


*/ 

*/ 

(sec)  */ 
*/ 
*/ 


/*  Flux  slope  for  Ga  ADJUST 

/*  Flux  slope  for  A13  ONLY  IF 

/*  Flux  slope  for  A14  CAL  DOES 

/*  Flux  slope  for  In  NOT  CONVERGE 

/*  Flux  slope  for  Sb  FAST  ENOUGH 

/*  Stabilization  time  for  Ga 

/*  Stabilization  time  for  A13 

/*  Stabilization  time  for  A14 

/*  Stabilization  time  for  In 

/*  Stabilization  time  for  Sb 


(prev  -11630)  */ 
*/ 

(prev  -13860)  */ 
(prev  -12080)  */ 


double 

TGamax 

_ 

999.  ; 

/*  Upper  Ga 

temp  limit 

*/ 

double 

TAl 3max 

= 

1249  . 

/* 

Upper 

A13 

temp 

limit 

*/ 

double 

TA14max 

- 

1249  . 

/* 

Upper 

A14 

temp 

limit 

*/ 

double 

T Inmax 

= 

899. 

/* 

Upper 

In 

temp 

limit 

*/ 

double 

TSbmax 

= 

799. 

/* 

Upper 

Sb 

temp 

limit 

*/ 

double 

TSimax 

= 

1399. 

/* 

Upper 

Si 

temp 

limit 

*1 

double 

TBemax 

= 

1149  . 

/* 

Upper 

Be 

temp 

limit 

*f 

double 

TGamin 

— 

600. 

/* 

Lower 

Ga 

temp 

limit 

*/ 

double 

TAl3min 

_ 

82  0. 

/* 

Lower 

A13 

temp 

limit 

*/ 

double 

TAl4min 

= 

820. 

/* 

Lower 

A14 

temp 

limit 

*/ 

double 

TInmin 

= 

400. 

/* 

Lower 

In 

temp 

limit 

*( 

double 

TSbmin 

= 

200. 

/* 

Lower 

Sb 

temp 

limit 

*/ 

double 

TSimin 

400. 

/* 

Lower 

Si 

temp 

limit 

*/ 

double 

TBemin 

= 

400. 

/* 

Lower 

Be 

temp 

limit 

*/ 

double  TGa_outg; 
double  TAl3_outg; 
double  TAl4_outg; 
double  TIn_outg; 
double  TSi_outg; 
double  TBe_outg; 
double  TSb_outg; 

double  fluxGa; 
double  f luxAl3 ; 
double  f luxAl4 ; 
double  fluxln; 
double  fluxSb; 


Outgas 

Outgas 

Outgas 

Outgas 

Outgas 

Outgas 

Outgas 


temp  for  Ga 
temp  for  A13 
temp  for  A14 
temp  for  In 
temp  for  Si 
temp  for  Be 
temp  for  Sb 


/*  Measured  flux  Ga 
/*  Measured  flux  A13 
/*  Measured  flux  Al4 
/*  Measured  flux  In 
/*  Measured  flux  Sb 


double  TGa; 

/* 

Temp  for  Ga 

*/ 

double  TAl 3 ; 

/* 

Temp  for  A13 

*/ 

double  TAl 4 ; 

/* 

Temp  for  A14 

.*/ 

double  Tin; 

/* 

Temp  for  In 

*/ 

double  TSb; 

/* 

Temp  for  Sb 

*/ 

double  errorGa; 

/* 

Ga  flux  error 

*/ 

double  errorAl3; 

/* 

A13  flux  error 

*/ 

double  errorAl4; 

/* 

A14  flux  error 

*/ 

double  errorln; 

/* 

In  flux  error 

*/ 

double  errorSb; 

/* 

Sb  flux  error 

*/ 

int  doneGa; 

/* 

Completion  flag 

for 

Gci 

(l=yes 

0  =  no) 

*/ 

int  doneAl3; 

/* 

Completion  flag 

for 

A13 

(l=yes 

0  =  no) 

*/ 

int  doneAl4; 

/* 

Completion  flag 

for 

A14 

(l=yes 

0  =  no) 

*/ 

int  doneln; 

/* 

Completion  flag 

for 

In 

(l=yes 

0  -  no) 

*/ 
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/*  Completion  £lag  £or  Sb  (l=yes  0  =  no) 


/ 


int  doneSb; 


int  Tbeamread  =  1; 
int  iread; 
int  attempt; 
int  log_idl; 


/*  Time  between  flux  readings  (sec)  */ 

/*  counter  during  flux  sampling  */ 

/*  Number  of  changes  of  cell  temp  */ 

/*  Data  logger  */ 


double  beam_flux; 
double  sum_flux; 
double  beam„f lux_open; 
double  beam_f lux_close ; 
double  background; 
long  tnow; 
long  tbegin; 
long  tleft;  . 

/**/ 

wait  until  start-up  time 


/*  Measured  flux  value  */ 
/*  Summation  variable  for  flux  calc  */ 
/*  Beam  flux  with  open  shutter  */ 
/*  Beam  flux  with  close  shutter  */ 
/*  Background  pressure  after  shutter  closed  */ 
/*  Present  time  */ 
/*  Time  to  start  up-ramp  */ 
/*  Time  left  before  start  of  up-ramp  */ 


/**/  ,  ' 
echo ( ) ; 

echo  {  "  =  =  :=  =  =  =  =  ““  =  =  =  =  =  :=  =  =  =  =  "  =  =  =  “  =  =  =:  =  =  =  =  =  =  =  ===  =  =  =  =  =  =  =  =  =  -  =  =  =  =  :j: 
echo ( "  l |  DON'T  FORGET  TO  PUT  SHUTTERS  AND  CAR  IN  REMOTE  | | 

echo  ( )  ; 


")  ; 
")  ; 
")  ; 


targ_time  =  hour  *  60  +  minute;  /*  Timer  code  from  J.  Vlcek  5/5/97  */ 

■fd  -  open(  "_clock",  0_RD0NLY  ); 


if  {  fd  <  0  )  { 

echo (  "I'm  sorry,  but  I'm  unable  to  access  the  system  clock  for  some  reason. 
echo(  "I  am  unable  to  schedule  your  command  file  for  later  execution."  ); 
echo (  "Please  send  an  email  to  software@epimbe.com  describing  this  problem, 
exit  {  EXIT__FAILURE  ); 

} 

echo  { )  ; 

echo (  "  The  up-ramp  is  now  scheduled  for  execution."  ); 


)  ; 


)  ; 


start_time  =  ioctl (  fd,  CLIOCTOD,  0  ) ; 

/*  Wait  until  midnight  if  the  target  time  is  earlier  in  the  day 
*  than  the  current  time  (ie,  the  file  executes  tomorrow) . 

*/ 

if  (  start_time  >=  targ_time  )  { 

while  (  ioctl (  fd,  CLIOCTOD,  0  )  >=  start_time  ) 

{ 

sleep {  20.0  ); 

> 

} 

while  (  ioctl {  fd,  CLIOCTOD,  0  )  <  targ_time  ) 

{ 

sleep (20.0); 

} 

close (  fd  ) ; 


/**/ 

/*==========  Test  chamber  pressure  before  cell  up  ramps 

/**/ 

echo{"  Testing  chamber  pressure"); 
load  ("pos3.cmd"); 

echo ( "  Wait  25  sec  for  substrate  to  turn  toward  cells*); 

echo ( ) ; 

sleep  (  25  ) ; 

iread  =  0; 

sum_f lux  =  0 ; 

while  (iread<10) 

{ 

beam_flux  -  reading { flux) ; 

echo ( "  Flux  gauge  =  ",  beam_flux) ; 

if(beam_flux  >  0) 

{ 

sum„flux  =  sum_flux  +  beam_flux; 
iread  =  iread+1; 

} 


=*/ 

=*/ 

=*/ 


/ *  Turn  CAR 


*/ 


Appendix 


sleep (  Tbeamread  ); 

} 

background  =  sum_f lux/10; 

if (  background  <  backtest  &&  background  >  1.0E-11  ) 

{ 

echo ( ) ; 

echo { "  The  chamber  pressure  ",  background  ); 

echo("  passed  the' test  limit  ",  backtest  ); 

echo ( "  The  cells  will  now  be  ramped  up  ") ; 
echo { ) ; 

} 

else 
{  , 

echo ( )  ; 

echo ( "  The  chamber  pressure  ",  background  ); 

echo{"  exceeds  the  test  limit  ",  backtest  ); 

echo{"  I  will  not  ramp  up  the  cells  "); 
echo ( ) ; 

kill {  getpid{ ) ,  SIGKILL) ; 

} 

/**/ 

/*===========  Ramp  up  cells  with  dummy  block  facing  cells  -----=======*/ 

/**/ 

/**/ 

/* -  Start  data  logger  - */ 

/**/ 

=  logger (20.0,  /*  log  every  20  seconds  */ 

’t ' , 

’ temp ( subs ) ' , 

' temp(Ga) ' , 

' temp (A13 ) '  , 

'temp (A14) ' , 

'temp (In) * , 

' temp ( Si ) ' , 

' temp (Be) ' , 

' is_open (Ga) ' , 

'is__open(Al3) ' , 

' is_open (A14) ' , 

'is_open(In) ' , 

' is_open (Si) ' , 

'is_open(Be) ' , 

'reading (flux) ' , 

"f luxcal . dat" ) ; 

/*  PUT  SB  BACK  WHEN  THE  EUROTHERM  IS  BACK,  */ 

/**/ 

/*— — _____  Turn  CAR  and  check  if  ramping  should  be  done 
/**/ 

load  ( "pos3 . cmd" ) ; 

echo("  Wait  25  sec  for  substrate  to  turn  toward  cells") 
echo ( ) ; 
sleep  (  25  )  ; 

if  (  rampGa  ==  1  j  |  rampAl3  =-  1  |  j  rampAl4  ==  1  | |  rarapln  - 
ramps i  ==  1  j  j  rampBe  ==  1  j  [  rampSb  ==  1  ) 

{ 

/**/ 

/* - -  Ramp  up  to  outgas  temperature  (2.5%  above  nominal  temp) 

/**/ 

echo ( "  Wait  ",  ramptime,  "  min  for  cells  to  ramp  up  "); 
echo  ( )  ? 

if  (  rampGa  ==  1  )  TGa  =  temp(  Ga  ); 

if  (  rampAl3  ==  1  )  TA13  =  temp{  A13  ); 

if  (  rampAl4  ==  1  )  TA14  =  temp(  A14  ); 

if  (  rampln  ==  1  )  Tin  =  temp (  In  ); 

if  (  rampSi  ==  1  )  TSi  =  temp(  Si  ); 

if  (  rampBe  ==  1  )  TBe  =  temp(  Be  ); 

if  (  rampSb  ==  1  )  TSb  -  temp (  Sb  ) ; 

TGa_outg  =  setpGa  *(1.0  +  outgas/100 . 0) ;  .  /*  Create  outgas  temp  */ 

TAl3_outg  =  setpA13* (1.0  +  outgas/100 . 0 ) ; 

TAl4_outg  =  setpAl4  * (1.0  +  outgas/100 . 0) ; 

TIn_outg  =  setpln  *(1.0  +  outgas/100 . 0) ; 

TSi_outg  =  setpSi  *(1.0  +  outgas/100 . 0) ; 


;=  1  |  j  /*  If  any  cell  is  to  be  ramped  */ 


/*  Get  current  setpoints  */ 


-*/ 

/*  Turn  CAR  */ 
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TBe_outg  =  setpBe  *(1.0  +  outgas/100 .0) ; 

TSb_outg  =  setpSb  *(1.0  +  outgas/100 .0) ; 

if  (  TGa_outg  >  TGamax  ) 

{ 

TGa„outg  =  TGamax; 

echo("  Warning  -  Ga  will  be  outgassed  at  max  temp  ",  TGamax,  "  C"); 

} 

if  (  TAl3_outg  >  TAl3max  ) 

{ 

TAl3_outg  =  TAl3max; 

echo("  Warning  -  A13  will  be  outgassed  at  max  temp  ",  TA13max, "  C")?' 
} 

if  (  TAl4_outg  >  TA14max  ) 

{ 

TAl4_outg  =  TA14max; 

echo("  Warning  -  A14  will  be  outgassed  at  max  temp  ",  TAl4max, "  C"); 

) 

if  {  TIn_outg  >  TInmax  ) 

{ 

TIn_outg  =  TInmax; 

echo("  Warning  -  In  will  be  outgassed  at  max  temp  ",  TInmax,"  C"); 

} 

if  (  TSi_outg  >  TSimax  ) 

{ 

TSi_outg  =  TSimax; 

echo("  Warning  -  Si  will  be  outgassed  at  max  temp  ",  TSimax,"  C"); 

} 

if  (  TBe_outg  >  TBemax  ) 

{TBe_outg  -  TBemax; 

echo { "  Warning  -  Be  will  be  outgassed  at  max  temp  ",  TBemax,"  C"); 

} 

if  (  TSb_outg  >  TSbmax  ) 

{ 

TSb_outg  =  TSbmax; 

echo("  Warning  -  Sb  will  be  outgassed  at  max  temp  ",  TSbmax,"  C " ) ; 

} 


if 

rampGa 

== 

1 

set_ramp ( 

Ga 

,  (TGa_outg 

-  TGa  ) /ramp time  ) ; 

/*  Set  new  ramp  rates 

*/ 

if 

rampAl3 

== 

1 

set_ramp ( 

A13 

, (TAl3_outg 

-  TA13) /ramptime  ); 

if 

rampAl4 

== 

1 

set_ramp ( 

A14 

,  (TAl4_outg 

-  TA14) /ramptime  ) ; 

if 

rampln 

=  = 

1 

set_ramp ( 

In 

,  (TIn_outg 

-Tin  ) /ramp time  ) ; 

/*  DEG /MIN 

*/ 

if 

rampSi 

=  = 

1 

set_ramp ( 

Si 

,  (TSi_outg 

-  TSi  ) /ramp time  ) ; 

if 

rampBe 

=  = 

1 

set_ramp ( 

Be 

, (TBe_outg 

-  TBe  ) /ramp time  ) ; 

if 

rampSb 

=  = 

1 

set_ramp ( 

Sb 

,  (TSb_outg 

-  TSb  ) /ramp time  ) ; 

if 

(  rampGa 

=  = 

1 

set_temp ( 

Gb. 

,TGa_outg) ; 

/*  Set  new  temperatures 

*/ 

if 

(  rampAl3 

== 

1 

set_temp ( 

A13 , TAl3_outg) 

if 

(  rampAl4 

=  = 

1 

set_temp( 

A14 

,TAl4_outg) 

if 

(  rampln 

=  = 

1 

set_temp( 

In 

, TIn_outg) ; 

if 

(  rampSi 

=  = 

1 

set_temp ( 

Si 

, TSi_outg) ; 

if 

(  rampBe 

=  = 

1 

set_temp ( 

Be 

,  TBe_outg) ; 

if 

(  rampSb 

=  = 

1 

set_temp ( 

Sb 

,  TSb_outg) ; 

sleep  (  ramptime*60) ;  /*  Wait  until  ramp  completed  */ 


/**/ 

/* -  Open  shutters  and  outgas  - */ 

/**/ 

echo("  Outgas  cells  " , outgastime,  "  min"}; 
echo  ( )  ; 

i f  (  rampGa  ==  1  )  shopen ( Ga  ) ; 

if  (  rampAl3  ==  1  )  shopen(Al3); 

if  (  rampAl4  ==  1  )  shopen (A14 ) ; 

if  (  rampln  ==  1  )  shopen (In  ) ; 

if  (  rampSi  ==  1  )  shopen (Si  ) ; 

if  (  rampBe  ==  1  )  shopen (Be  ); 

if  (  rampSb  ==  1  )  shopen (Sb  ); 

sleep  {  outgastime*60) ; 

/**/ 

/* -  Close  shutters  and  ramp  down  to  setpoints  - */ 

/**/ 

echo("  Close  cells  and  ramp  to  setpoints  wait  5  min"); 
echo ( ) ; 

if  (  rampGa  ==  1  }  shclose(Ga  );  /*  Close  the  shutters  */ 

if  (  rampAl3  ==  1  )  shclose (A13 ) ; 

if  (  rampAl4  ==  1  )  shclose (A14) ; 
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if  (  rampln  ==  1  )  shclosedn  ); 

if  (  rampSi  ==  1  )  shclose(Si  ); 

if  (  rampBe  =-  1  )  shclose(Be  ); 

if  (  rampSb  ==  1  )  shclose(Sb  ); 

if  {  rampGa  ==  1  )  TGa  =  temp (  Ga  };  /*  Get  current  setpoints  */ 

if  (  rampAl.3  ==  1  )  TA13  =  temp  (  A13  ); 

if  {  rampA!4  ==  1  )  TA14  =  temp(  A14  ); 


if  i 

f  rampln 

==  i  ; 

I  Tin 

=  tempf 

In 

if  < 

(  rampSi 

=-  l  ; 

1  TSi 

=  tempf 

Si 

if  i 

(  rampBe 

==  l  ; 

\  TBe 

-  tempf 

Be 

if  < 

(  rampSb 

==  l  ] 

1  TSb 

=  tempf 

Sb 

if  (  rampGa  ==  1  )  set„ramp (  Ga  , (TGa  -  setpGa  )/5.  );  /*  Set  new  ramp  rates  */ 

if  {  rampAl3  ==  1  )  set_ramp (  A13,(TA13  -  setpAl3)/5.  ); 

if  (  rampAl4  ==  1  )  set_ramp (  A14, (TA14  -  setpA14)/5.  ); 

if  (  rampln  =-  1  )  set_ramp (  In  , (Tin  -  setpln  )/5.  )  ;.  /*  DEG /MIN  */ 

if  (  rampSi  ==  1  )  set_ramp (  Si  , {TSi  -  setpSi  )/5.  )? 

if  (  rampBe  ==  1  )  set_ramp (  Be  , {TBe  -  setpBe  )/5.  ); 

if  {  rampSb  ==  1  )  set_ramp (  Sb  , (TSb  -  setpSb  )/5.  ); 

if  (  rampGa  ==  1  )  set_temp (  Ga  , setpGa  );  /*  Set  new  temperatures  */ 

if  {  rampAl3  ==  1  )  set_temp(  Al3,setpAl3  ) ; 

if  {  rampAl4  ==  1  )  set_temp (  Al4,setpAl4  ) ; 

if  (  rampln  ==  1  )  set_temp (  In  , setpln  ); 

if  (  rampSi  -=  1  )  set_temp (  Si  , setpSi  ); 

if  (  rampBe  ==  1  )  set_temp (  Be  , setpBe  ); 

if  (  rampSb  ==  1  )  set_temp (  Sb  , setpSb  ); 

sleep  (  300  } ; 

/**/ 

/* -  Wait  5  more  minutes  to  ensure  stability  - */ 

/**/ 

echo("  Wait  5  min  for  stability  "); 
echo { ) ; 

sleep  (  300  ) ; 

}  /*  End  of  cell  excercise  */ 

/**/ 

/*===========;=====  Measure  Ga  flux  =========================*/ 

/**/ 

iff  calGa  ==  1  ) 

{ 

doneGa  =  0; 
attempt  =  0; 

echof”  Measure  Ga  flux"); 
echo ( ) ; 

TGa  =  setpGa; 

set_ramp (Ga, 100) ;  /*  set  fast  rate  for  small  adjustments  */ 

} 

while  {  doneGa  ==  0  &&  calGa  ==  1  ) 

{ 

1**1 

f* -  Set  new  temperature  turn  flux  guage  away  from  cells  and  wait  for  stabilization  - */ 

/**/ 

iff  TGa  <  TGamax  &&  TGa  >  TGamin  ) 
set_temp ( Ga , TGa ) ; 
else 
{ 

echof"  Ga  setpoint  outside  allowd  interval  -  Process  terminated"); 
kill (  getpidf),  SIGKILL) ; 

} 

attempt  =  attempt  +  1; 

echof"  Seting  new  temp  =  ",TGa,"  and  waiting  ",stabGa,"  sec."); 
echof"  Time  is:  ",mctime (time (0) )  ); 
echo  f ) ; 

load  ( "pos3 . cmd" ) ; 
sleep  {  stabGa  ) ; 

/**/ 

/* -  ,  open  shutter  and  wait  — - ■•*/  . 

1**1 

shopen(Ga) ; 

echof"  Wait  ",  transient , "  sec  during  transient"); 
echo ( ) ? 

sleep  (  transient  ) ; 
load  { "posl . cmd" ) ; 
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echo("  Wait  30  sec  for  guage  to  turn  toward  cells"); 
echo  ( )  ; 
sleep  {  30  )  ; 

/**/ 

/* -  Measure  with  shutter  open  - */ 

/**/ 

echo("  Measure  with  shutter  open")? 

echo ( ) ; 

iread  =  0; 

sum_flux  -  0; 

while  (iread<10) 

{ 

beam_flux  =  reading ( flux) ; 
echo("  Flux  gauge  =  ",  beam_flux) ? 
if (beam_f lux  >  0) 

{ 

suin_flux  -  sum_f lux  +  beam_f lux; 
iread  =  iread+1; 

} 

sleep (  Tbeararead  )? 

} 

beam_f lux_open  =  sum_flux/10; 
echo ( ) ; 

echo{"  Average  Flux  =  ",  beara_flux_open) ; 
echo { ) ; 

/**/ 

/* -  Measure  with  shutter  closed  - */ 

/**/ 

shclose (Ga) ; 
sleep (20) ; 
iread  =  0; 
sum_flux  -  0; 
while  (iread<10) 

{ 

back_flux  =  reading (flux) ; 

echo ( "  Flux  gauge  =  ",  back_flux) ? 

if(back_flux  >  0  ) 

{ 

sum_flux  =  sum_flux  +  back_flux; 
iread  =  iread+1? 

} 

sleep (  Tbeamread  ) ? 

) 

beam_f lux_close  =  sum_f lux/10 ; 
echo  { )  ; 

echo("  Average  Flux  =  ",  beam_flux_close) ? 

■  echo ( ) ; 

fluxGa  =  beam_flux_open  -  beam_f lux_close; 
echo{"  Net  Flux  =  ",  fluxGa); 
echo  ( )  ? 

/**/ 

I* -  Test  flux  and  calc  temp  correction - */ 

/**/ 

errorGa  =  (  f luxGa-targetGa  )/targetGa; 
echo ( "  Ga  error  =  " , errorGa  ) ? 
if (  fabs{  errorGa  )  >  precGa  ) 

{ 

TGa  =  1. / (1 . / (TGa+273)  -  ( loglO ( fluxGa) -loglO ( targetGa )) /slopeGa  )  -  273; 

) 

else 

doneGa  =  1 ; 

} 

/**/ 

/* -  Ga  calibrated  - */ 

/**/ 

if(  calGa  ==  1  ) 

^set  ramp(Ga,10);  /*  reset  slow  rate  for  protection  */ 

echo ( ) ; 

echo("  Ga  calibration  converged  in  ", attempt,"  attempts")? 
echo("  Final  error  was  ", errorGa  ); 
echo ( ) ; 

) 

/**/ 
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/*z=====s===5S====  End  0f  Ga  loop  =:“======================:*/ 

/**/ 

/**/ 

/*==============-=  Measure  A13  flux  ================—=======*/ 

/**/ 

iff  calAl3  ==  1  ) 

{ 

doneAl3  =0; 
attempt  =  0; 

echof”  Measure  A13  flux"); 
echo ( ) ; 

TA13  =  setpAl3; 

set_ramp (A13, 100) ;  /*  set  fast  rate  for  small  adjustments  */ 

} 

while  (  doneAl3  ==  0  &&  calAl3  ==  1  ) 

{ 

/**/ 

/* -  Set  new  temperature,  turn  flux  guage  away  from  cells  and  wait  for  stabilization  - */ 

/**/ 

iff  TA13  <  TAl3max  &&  TA13  >  TA13min  ) 
set_temp (A13 , TA13 ) ; 
else 
{ 

echof”  A13  setpoint  outside  allowd  interval  -  Process  terminated"); 
kill (  getpidf),  SIGKILL) ; 

} 

attempt  =  attempt  +  1; 

echof”  Seting  new  temp  =  ”,TA13,"  and  waiting  ”,stabAl3,"  sec"); 
echof"  Time  is:  " ,mctime (time (0) )  ); 
echo ( ) ; 

load  ( "pos3 .cmd" ) ; 
sleep  {  stabAl3  ) ; 

/**/ 

/* -  ,  open  shutter  and  wait  — - */ 

/**/ 

shopen  (A13 )  ; 

echof”  Wait  ”,  transient , "  sec  during  transient"); 
echo  ( )  ; 

sleep  (  transient  ) ; 
load  ( ”posl . cmd" ) ; 

echof”  Wait  30  sec  for  guage  to  turn  toward  cells"); 
echo  { )  ; 
sleep  (  30  )  ; 


/* -  Measure  with  shutter  open  - */ 

/**/ 

echof”  Measure  with  shutter  open"); 

echo ( ) ; 

iread  =  0; 

sum_f lux  =  0; 

while  (ireaddO) 

{ 

beam„flux  =  reading ( flux) ; 
echof”  Flux  gauge  =  ”,  beam_f lux) ; 
if (beam_flux  >  0) 

{ 

sum_f lux  =  sum_flux  +  beam_flux; 
iread  =  iread+1; 

} 

sleep (  Tbeamread  ); 

} 

beam_f lux_open  -  sum_f lux/10; 
echo ( ) ; 

echof”  Average  Flux  =  ”,  beam_flux_open) ; 
echo  ( )  ; 

/**/ 

/* -  Measure  with  shutter  closed  - */ 

/**/ 

shclose (A13) ; 
sleep (20) ; 
iread  =  0; 
sum_f lux  =  0; 
while  (iread<10) 
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{ 

back_flux  =  reading ( flux) ; 
echo{"  Flux  gauge  =  ",  back_f lux) ; 
if(back_flux  >  0  ) 

{ 

sum„flux  -  sum„flux  +  back_flux,- 
iread  =  iread+1; 

} 

sleep (  Tbeamread  ); 

} 

beam_f lux_close  =  sum_f lux/10; 
echo  U  ; 

echo{"  Average  Flux  =  ",  beam_flux_close) ; 
echo  ( )  ; 

fluxAl3  ss  beam_f lux_open  -  beam_f lux_close; 
echo("  Net  Flux  =  ",  fluxAl3) ; 
echo ( ) ; 

/**/ 

/* -  Test  flux  and  calc  temp  correction - */ 

/**/ 

errorAl3  =  (  f luxAl3-targetAl3  )/targetAl3; 
echo("  A13  error  =  ",errorAl3  ); 
if(  f abs  (  errorAl3  .)  >  precAl3  ) 

{ 

TA13  =  1 . / (1 . / (TA13+273)  -  { loglO { f luxAl3 ) -loglO ( targetAl3 ) ) /slopeAl3  )  -  273; 

} 

else 

doneAl3  =  1; 

} 

/**/ 

/* -  AX 3  calibrated  - */ 

/**/ 

if (  calAl3  ==  1  ) 

{ 

set_ramp (A13 , 10) ;  /*  reset  slow  rate  for  protection  */ 

echo { ) ; 

echo("  A13  calibration  converged  in  ", attempt,"  attempts"); 
echo{"  Final  error  was  ",errorAl3  ); 
echo ( ) ; 

} 

/**/ 

/*;=======:=:========  End  of  Al3  loop  ========*““=======““=*/ 

/**/ 

/*================  Measure  A14  flux  ===========■==============* / 

/**/ 

if  (  calAl4  ==  1  ) 

{ 

doneAl4  =  0; 
attempt  =  0; 

echo("  Measure  A14  flux"); 
echo ( ) ; 

TA14  =  setpAl4 ; 

set__ramp{Al4, 100) ;  /*  set  fast  rate  for  small  adjustments  */ 

} 

while  (  doneAl4  ==  0  &&  calA14  ==  1  ) 

{ 

/**/ 

/* -  Set  new  temperature,  turn  flux  guage  away  from  cells  and  wait  for  stabilization  - */ 

/**/ 

if (  TA14  <  TAl4max  &&  TA14  >  TAl4min  ) 
set_temp (A14 , TA14 ) ; 
else 
{ 

echo("  A14  setpoint  outside  allowd  interval  -  Process  terminated"); 
kill  (  getpidO,  SIGKILL)  ; 

} 

attempt  =  attempt  +  1; 

echo{"  Seting  new  temp  =  ",TA14,"  and  waiting  ",stabAl4,"  sec"}; 
echo{"  Time  is:  " ,mctime (time (0) )  ); 
echo  ( )  ; 

load  ("pos3.cmd"); 
sleep  {  stabA14  ) ; 
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1**1 

/* - -  open  shutter  and  wait  - * / 

/**/ 

shopen(A14) ; 

echo ( "  Wait  ",  transient , "  sec  during  transient"); 
echo ( ) ; 

sleep  (  transient  )  ; 
load  { "posl .cmd" ) ; 

echo("  Wait  30  sec  for  guage  to  turn  toward  cells"); 
echo ( ) ; 
sleep  (  30  ) ; 

/**/ 

/* -  Measure  with  shutter  open  - 

/**/ 

echo { "  Measure  with  shutter  open"); 
echo ( ) ; 
iread  =  0; 
sum_f lux  =  0; 
while  (iread<10) 

{ 

beam_flux  =  reading ( flux) ; 
echo{"  Flux  gauge  =  ",  beam_flux) ; 
if (beam_f lux  >  0) 

{ 

sum_flux  =  sum_flux  +  beam_flux; 
iread  =  i'read+1 ; 

> 

sleep (  Tbeamread  ); 

} 

beam_f lux_open  -  sum_flux/10; 
echo ( ) ; 

echo ( "Average  Flux  =  ",  beam_flux_open) ; 
echo { ) ; 

1**1 

f* -  Measure  with  shutter  closed  - 

1**1 

shclose (A14 ) ; 
sleep (20) ; 
iread  =  0; 
sum_f lux  =  0; 
while  (iread<10) 

{ 

back_flux  =  reading ( flux) ; 
echo{"  Flux  gauge  -  ",  back_f lux) ; 
if (back_f lux  >  0  ) 

{ 

sum_flux  =  sum_flux  +  back_flux; 
iread  =  iread+1; 

} 

sleep (  Tbeamread  ) ; 

} 

beam_f lux_close  =  sum_f lux/10; 
echo { ) ; 

echo("  Average  Flux  =  ",  beam_flux_close) ; 
echo ( ) ; 

fluxAl4  =  beam_f lux_open  -  beam_f lux_close; 
echo("  Net  Flux  =  ",  fluxAl4) ; 
echo { ) ; 

1**1 

/* -  Test  flux  and  calc  temp  correction - */ 

/**/ 

errorAl4  =  (  f luxAl4-targetAl4  )/targetAl4; 
echo{"  A14  error  =  ",errorAl4  ); 
if (  f abs (  errorAl4  )  >  precAl4  ) 

{ 

TA14  =  1 . / ( 1 . / (TA14+273 )  -  ( loglO { fluxAl4 ) -loglO { targetAl4 ) ) /slopeAl4  )  -  273 

} 

else 

doneA14  =  1; 

} 

/**/ 

/* -  A14  calibrated  - */ 

1**1 

if {  calA14  ==  1  ) 

{ 


■*/ 


■*/ 
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set_ramp(A14, 10) ;  /*  reset  slow  rate  for  protection  */ 

echo ( )  ; 

echo{”  A14  calibration  converged  in  ”,  attempt,"  attempts"); 
echo(”  Final  error  was  ”,errorAl4  ); 
echo  ( )  ; 

} 

/**/ 

End  of  A14  loop  ========:“=======s========-*/ 

/**/  *t 
/*=====:==:=========  Measure  In  flux  ===s=====================:*  / 

/**/ 

if (  calln  ==  1  )  > 

{ 

doneln  -  0; 
attempt  =  0; 

echo("  Measure  In  flux")? 
echo ( ) ; 

Tin  =  setpln; 

Set_ramp ( In, 100 ) ;  •  /*  set  fast  rate  for  small  adjustments  */ 

) 

while  {  doneln  ==  0  &&  calln  --  1  ) 

{ 

/**/ 

/* -  Set  new  temperature,  turn  flux  guage  away  from  cells  and  wait  for  stabilization  */ 

/**/ 

if(  Tin  <  TInmax  &&  Tin  >  TInmin  ) 
set_temp (In, Tin) ; 
else 
{ 

echo(”  In  setpoint  outside  allowd  interval  -  Process  terminated")? 
kill  (  getpidO,  SIGKILL)  ; 

} 

attempt  =  attempt  +  1; 

echo ( ”  Seting  new  temp  =  ”,TIn,"  and  waiting  ”,stabln,"  sec"); 
echo(”  Time  is:  ” ,mctime (time (0) )  ); 
echo { ) ; 

load  ( ”pos3 -cmd" ) ; 
sleep  (  stabln  ) ; 

1**1 

/* -  open  shutter  and  wait  - */ 

I**/ 

shopen(In) ; 

echo { ”  Wait  ”,  transient , "  sec  during  transient"); 
echo ( ) ; 

sleep  (  transient  ) ; 
load  (”posl .cmd" ) ; 

echo(”  Wait  30  sec  for  guage  to  turn  toward  cells"); 
echo ( ) ; 
sleep  (  30  ) ; 

I**/ 

/* -  Measure  with  shutter  open  - */ 

/**/ 

echo(”  Measure  with  shutter  open"); 

echo  ( )  ; 

iread  =  0; 

sum_flux  “  0; 

while  (iread<10) 

{ 

beam_flux  -  reading (flux) ; 
echo(”  Flux  gauge  =  ”,  beam_f lux) ; 
if (beam_f lux  >  0) 

{ 

sum_flux  =  sum_flux  +  beam_flux? 
iread  =  iread+1; 

} 

sleep (  Tbeamread  ) ; 

} 

beam_f lux_open  =  sum_f lux/10; 
echo  ( )  ; 

echo { ”  Average  Flux  =  ”,  beam_flux_open) ; 
echo { )  ; 

/**/ 

/* -  Measure  with  shutter  closed  - */ 
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shclose (In) ; 
sleep (20) ; 
iread  =  0; 
sum_f lux  =  0; 
while  (iread<10) 

{ 

back_flux  =  reading ( flux) ; 
echo("  Flux  gauge  =  ”,  back_f lux) ; 
if(back_flux  >  0  ) 

{ 

sum__flux  =  sum_flux  +  back_flux; 
iread  =  iread+1; 

) 

sleep (  Tbeamread  ); 

} 

beam_f lux„close  =  sum_f lux/10; 
echo ( ) ; 

echo(”  Average  Flux  =  ",  beam_flux_close) ; 
echo ( ) ; 


fluxln  =  beam_f lux_open  -  beam_f lux_close; 
echo("  Net  Flux  =  ",  fluxln); 
echo  ( )  ; 

/**/ 

/* -  Test  flux  and  calc  temp  correction - - - */ 

/**/ 


errorln  =  (  f luxln-targetln  )/targetln; 
echo ( "  In  error  =  " , errorln  ) ; 
if(  f abs (  errorln  )  >  precln  ) 

{ 

Tin  =  l./(l./(TIn+273)  -  (loglO ( fluxln) -loglO (targetln) } /slopeln  )  -  273; 

> 

else 

doneln  =  1; 

} 

/**/ 

/* -  in  calibrated  - */ 


/**/ 

if(  calln  ==  1  ) 

{ 

set_rarap ( In, 10 ) ; 
echo  ( )  ; 

echo ("In  calibration  converged  in  ",  attempt,"  attempts"); 
echo ("Final  error  was  ”, errorln  ); 
echo ( ) ; 

} 

/**/ 

/*======:=  =  =  =  === ====  =  ===  ===  ======  =  ===  =  ====  =  =  =  ===  ====  =  =  =  =  =  =  =  =  ===  = 

En(j  q£  jn  3_00p  ==  =  ===  ===  =  =  =  =  =  =  === 

/*============================================================ 

/**/ 

/*============================================================ 


f*  reset  slow  rate  for  protection 


*/ 

:*/ 

:*/ 

*/ 


/*================  Measure  Sb  flux 

/*=====================================: 


*/ 

*/ 


/**/ 


*/ 


if (  calSb  ==  1  ) 

'{ 

doneSb  =  0; 
attempt  -  0; 

echo ( ”  Measure  Sb  flux"); 
echo ( ) ; 

TSb  =  setpSb; 

set_ramp (Sb, 100) ;  /*  set  fast  rate  for  small  adjustments 

} 

while  (  doneSb  ==  0  &&  calSb  ==  1  ) 

{ 


/**/ 


/* -  Set  new  temperature,  turn  flux  guage  away  from  cells  and  wait  for  stabilization  - */ 

/**/ 

if (  TSb  <  TSbmax  &&  TSb  >  TSbmin  ) 
set_temp (Sb, TSb) ; 
else 


{ 

echo{"  Sb  setpoint  outside  allowd  interval  -  Process  terminated"); 
kill  (  getpidO,  SIGKILL)  ; 

} 

attempt  =  attempt  +  1; 


*/ 
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echo ( "  Seting  new  temp  =  ",TSb, "  and  waiting  ",stabSb,"  sec"); 
echo("  Time  is:  " ,mc time (time (0) )  ); 
echo  ( )  ; 

load  ( "pos3 . cmd" ) ; 
sleep  (  stabSb  ) ; 

/**/ 

/* -  open  shutter  and  wait  - */ 

1**1 

shopen(Sb) ; 

echo("  Wait  ",  transient , "  sec  during  transient"); 
echo ( ) ; 

sleep  (  transient  ) ; 
load  ( "posl . cmd" ) ; 

echo("  Wait  30  sec  for  guage  to  turn  toward  cells"); 
echo ( ) ; 
sleep  {  30  )  ; 

/**/ 

/* -  Measure  with  shutter  open  - - */ 

/**/ 

echo { *  Measure  with  shutter  open"); 

echo  ( )  ; 

iread  -  0; 

sum_flux  =  0; 

while  (iread<10) 

{ 

beam_flux  =  reading ( flux) ; 
echo("  Flux  gauge  =  ",  beam_f lux) ; 
if(beam_flux  >  0) 

{ 

sum_flux  =  sum_flux  +  beam_flux; 
iread  =  iread+1; 

} 

sleep (  Tbeamread  ); 

} 

beam_f lux_open  =  sum_flux/10; 
echo { ) ; 

echo{"  Average  Flux  =  ",  beam_f lux_open) ; 
echo ( ) ; 

/**/ 

/* -  Measure  with  shutter  closed  - */ 

/**/ 

she lose (Sb) ; 
sleep (20) ; 
iread  =  0; 
sum_flux  -  0; 
while  (iread<10) 

{ 

back_flux  =  reading ( flux) ; 
echo("  Flux  gauge  =  ",  back_flux) ; 
if(back_flux  >  0  ) 

{ 

sum_flux  =  sum_flux  +  back_flux; 
iread  =  iread+1; 

} 

sleep (  Tbeamread  ); 

) 

beam_f lux_close  =  sum_f lux/10 ; 
echo ( ) ; 

echo("  Average  Flux  =  ",  beam_f lux_close) ; 
echo ( ) ; 

fluxSb  =  beam_f lux_open  -  beam_f lux_close; 
echo("  Net  Flux  =  ",  fluxSb); 
echo ( ) ; 

/**/ 

/* -  Test  flux  and  calc  temp  correction - */ 

/**/ 

errorSb  =  (  f luxSb-targetSb  )/targetSb; 
echo ( "  Sb  error  =  " , errorSb  ) ; 
if (  f abs (  errorSb  )  >  precSb  ) 

{ 

TSb  =  1 . / (1 . / (TSb+273)  -  ( loglO ( fluxSb) -loglO { targetSb) ) /slopeSb  ) 
} 

else 

doneSb  =  1; 


-  273; 
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End  of  Sb  loop 


/**/ 

/* -  Sb  calibrated  - */ 

/**/ 

if(  calSb  ==  1  ) 

{ 

set_ramp {Sb, 10); 
echo ( ) ; 

echof"  Sb  calibration  converged  in  ",  attempt,"  attempts"); 
echo ( "  Final  error  was  "  ,  errorSb  ) ; 
echo { ) ; 

} 

/**/ 

/*============================================================= 

/*================ 

/*============================== 

/**/ 

/**/ 

/*====================^=s=^=== 

/*==-=============  Summary 

/**/ 

kill (log_idl,  SIGTERM) ; 
echo { ) ; 

if (  calGa  ~  1  }  echof" 
iff  calAl3  ==  1  )  echof" 
iff  calAl4  ==  1  )  echof" 
iff  calln  ==  1  )  echof" 
iff  calSb  ==  1  )  echof" 
echo  f ) ; 

echof"  A  record  of  the  temperatures,  shutter  status  and  flux  values 
echof"  is  stored  in  the  file  FLUXCAL.DAT  "); 
echo ( ) ; 


/*  Stop  logging  */ 


/*  reset  slow  rate  for  protection  */ 


Final 

Ga 

error 

was 

* , errorGa  , 

"  TGa  =  " , TGa  ) ? 

Final 

A13 

error 

was 

" , errorAl3 , 

"  TA13  =  " , TA13 ) ; 

Final 

A14 

error 

was 

" , errorAl4 , 

"  TA14  =  " , TAl 4 ) ; 

Final 

In 

error 

was 

" , errorln  , 

"  Tin  =  " , Tin  ) ; 

Final 

Sb 

error 

was 

" , errorSb  , 

"  TSb  =  " , TSb  ) ; 
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